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ABSTRACT 

It is demonstrated that improved resolution measurements of pre- 
cipitation may be obtained from satellite borne radars with antenna 
beams having relatively large surface footprints. The method em- 
ploys deconvolution and Fourier transform procedures, and assumes a 
knowledge of the antenna beam pattern. As an example, the technique 
is specifically directed towards the application of future space- 
borne radar altimeters which may contain additional range gates to 
enable the measurement of rain at altitude. It is demonstrated that 
because of the natural variability of rain in the lateral extent, 
the standard beam averaging over the footprint could easily produce 
erroneous interpretations of the intensity of rain and its extent. 

On the other hand, many of these ambiguities may be removed employing 
the deconvolution techniques described. 

Rain measurement of the type described here are considered 
vital from the standpoint of representing a flag for altimeter data 
that may be corrupted by rain. It also provides sorely needed rain 
data over the oceans where little or no such data is available. 

1. INTRODUCTION 

It has been pointed out that rain may distort and/or attenuate 
the sea surface echo returns of signals originating from spaceborne 
radar altimeters, resulting in gross errors of the mean sea height 
(Walsh, et al., 1983). The employment of a "piggyback" modification 
to future spaceborne altimeter systems has been suggested for mea- 
suring or monitoring the presence of rain at altitude (Goldhirsh and 
Walsh, 1982) . The experimenter may employ this rain information as 
a flag for culling out the concomitant ocean surface measurements, 
as well as to provide additional meteorological information regard- 
ing the characteristics of precipitation. 

A difficulty that is encountered in the measurement of precip- 
itation with satellite borne, range gated radars deals with the poor 
resolution encountered because of the large beamwidth associated 
with the radar antenna located at relatively large distances from 
the earth. For example, the 3 dB down points for the one way an- 
tenna gain function corresponding to the Seasat altimeter antenna 
(1.6 degree beamwidth) located at an altitude of 800 km, gave rise 
to a 22.3 km diameter footprint at the earth's surface. Rain is, 
in general, highly variable along the horizontal extent, especially 
for convective types of precipitation (e.g., thunderstorms). Even 
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for stratiform rains (e.g., widespread precipitation), it is common 
for imbedded or convective cells of more intense rain rates to 
exist (Goldhirsh, 1983) . As the return echo signal from a precip- 
itating medium represents the convolution of the square of the gain 
and the true rain reflectivity profile (as will be demonstrated in 
Section 2) , the rain reflectivity deduced from the beam averaged 
echo signal may be significantly smaller than the true reflectivity 
level. Furthermore, as the antenna sweeps by the rain structure on 
the ground, the inevitable beam averaging may smooth out the highly 
variable rain structure such as to give a distorted picture of the 
ground rain intensity profile. 


2. FORMULATION OF THE PROBLEM 


Consider the configuration depicting a satellite at an altitude, 
Rq, above the mean sea surface. The satellite is assumed to have a 
pulsed radar having a downward pointing antenna with a beamwidth, 0 O . 
We shall examine here the backscatter from the precipitation con- 
tained within a pulse volume centered at a range R. 

The echo power is primarily due to backscatter from an array of 
raindrops within a pulse volume defined by the beamwidth and the 
range resolution of the pulse. 

This echo power may be shown to be given by (Goldhirsh and 
Monaldo, 1983), 


P = C • Z 
roc 


where 


Z (x ) 
c o 


where , 


C c - [( liillni ) V Vl F1B) V-R ^77] 


dxdy 


( 1 ) 

( 2 ) 

(3) 


o o 

where |k 0 | = | (nr -1) / (in +2) | and where the quantity Z c (x Q ) is refer- 

red to as the effective beam averaged radar reflectivity factor where 
C Q is an effective constant. Also, m is the complex reflective index 
of water drops, c is the velocity of light (m/sec) , T is the pulse 
width (sec), Pt is the transmitted power (watts), X is the wavelength 
(meters), Lt*L r is the transmitter and receiver losses, respectively 
(£1) , R is the range from the radar to the center of the pulse vol- 
ume, x^,y^ are the footprint lengths along the principal planes of 
the gain function, x Q is the displacement of the beam nadir along the 
x direction in units of meters from a fixed origin at the ground and 
x,y is the planar dimensions relative to a fixed origin. 


In the derivation of (2) we have integrated over the range res- 
olution in the range direction given by 6 = (cT/2)*F(B), where F(B) 
is an effective loss factor due to an assumed Gaussian filter re- 
sponse at the receiver (Doviak and Zrnic, 1979). For a matched fil- 
ter linear receiver, F (B) * -2.3 dB. As a representative beam struc- 
ture, we have assumed in the derivation of (2) the Probert-Jones 
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(1962) form of the one way gain function; namely a Gaussian struc- 
ture. 


Z €a (x,y) is the effective radar reflectivity factor for rain 
defined by, 


Z 


eq 



0(D) N (D) dD 


(m 3 ) 


( 4 ) 


where, o(D) is the backscatter cross section for a drop of diameter, 
D and N(D)dD is the number of drops between the sizes of D and D^D 
per unit volume. 


It is also tacitly assumed that the attenuating effects due to 
precipitation may be ignored. This assumption becomes tenable if the 
pulse volume near the top of the rain region is sampled. This may 
create the condition such that there is a negligible path length over 
which the signal propagates within the sampled volume. 


3. BEAM AVERAGED RADAR REFLECTIVITY FACTOR EXPRESSED AS A CONVOLU- 
TION INTEGRAL 

Assume the equivalent reflectivity factor within the integrand 
of (2) may be expressed by the product relationship 


Z eq = Z l (x) Z 2 (y) 
Substituting (5) into (2) 


where 


where 


/ 81n2 \ 
‘ \ x l y l/ 


I I (x ) 
Y o 


uo 

I„ = / Z 2 (y) g (y) dy 

— oo 

00 

I(x q ) = J Zj^x) g(x-x Q ) dx 


g(y) = exp 
g(x) * exp 


- (81n2) 


- (81n2) 


ft)*] 

(kf\ 


(5) 


( 6 ) 

(7) 

( 8 ) 


(9) 

( 10 ) 


It is interesting to note that I (x ) as given by (8) may alter- 
nately be expressed as the convolution of the gain square function, 
g (x) , and the reflectivity factor, Z^x), where both are taken in the 
x direction. Hence, 


I (x Q ) « Zj^x) * g (x) 


(ID 
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where the * denotes the convolution. 

4. EXTRACTING THE RADAR REFLECTIVITY FACTOR PROFILE FROM THE CONVO- 
LUTION INTEGRAL 

By taking the Fourier transform of both sides of (11) , and em- 
ploying a fundamental property of the convolution integral , we obtain 


where 


I(m) 

Z 1 (u) 

G(w) 



= Z^(oj) G(w) 

I (x Q ) exp(-ju)X Q ) dx Q 
Z^(x) exp(-jux) dx 




exp(-jujx) dx 


(12) 

(13) 

(14) 

(15) 


where I (to), Z^(to), and G(to) are respectively the Fourier transforms 
of I (x Q ) , Zi(x) and g (x) are given by (13) through (15). It is ap- 
parent that by solving for Z^((o) in (12) and applying the inverse 
transform to both sides (i.e., FT“1 {}), the profile Z. (x ) may be 
established. That is, 0 

00 

W - FT ' 1 {§{SfJ - 57 / d “ < 16 > 

— 00 


Equation (16) represents the general formula for extracting the 
reflectivity profile in the x direction, given a knowledge of the 
measured echo power, P r (x c ), obtained as the beam sweeps the rain 
profile as a function of x Q . 

Substituting (10) into (15) , the Fourier transform of the gain 
function is given by, 


1/2 


G(w) » ± 


ITXj 

rnr? 


2 2 
h> X^ 


exp I" OTnJ 

Substituting (17) into (16), we obtain, 


(17) 


2 l (x o> 


00 

/ 2* *0 f 

CM 

X 

CM 

3 

pT* J lM “ p 

[ + mnl + jwx oJ 

'* *1 ’ — 



du 


(18) 


Equation (18' represents the resulting expression for the reflec- 
tivity profile, where I (w) is the Fourier transform of I(x ) obtained 
through the beam averaged power measurement. 
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5. EXPECTED POWER LEVELS FOR A REPRESENTATIVE SATELLITE BORNE RADAR 
General Case 


In this Section we estimate the minimum detectable beam aver- 
aged reflectivity factor and the corresponding rain rate. As an 
example, we shall assume similar radar parameters as those for the 
Seasat radar altimeter (MacArthur, 1978) . These are given as fol- 
lows: A = 2.22 x 10“2 m; P t = 2 x 10-3 watts; T = 3.2 x 10”6 sec; 

R = 8 x 10 5 m; 0 L = = 1.6° (2.792 x 10“ 2 radians); G = 40.6 dB 

(1.148 x 10 4 ); L T = -0.9 dB (0.813); Lr = -1.2 dB (0.758); | K 0 | 2 = 
0.9; F(B) = -2.3 dB (0.59). Substituting these parameter values 
into (3) , the power reduces to 

P r (dBm) = - 140.5 + dBZ (19) 

6 3 

where the reflectivity factor, Z c is in mm /m , P is in dBm, and 
C = 8.9 x 10~15. Also dBZ is defined by, 

dBZ = 10 Log ., q Z c (20) 

Assuming a unity signal to noise ratio of -115 dBm (5 dB smaller than 
the Seasat case), and solving for dBZ, we obtain (dBZ) m ^ n = -115 + 
140.5 = 25.5. Hence, the radar is capable of detecting beam averaged 
reflectivity factors of 25.5 dB and larger with a signal to noise 
ratio of 1 or larger. 


To further amplify the above example, we assume the following 
empirical relation relating the rain rate P expressed in mm/hr, 
to the reflectivity factor, expressed in mm^/m^, 

R r " <Z c /200) 0,625 (21) 


This relationship is consistent with the assumption of a Marshall- 
Palmer drop size distribution (1948) . Relating (21) and (20) 


R r = (3.65 x 10~ 2 ) 10 <* 0625 ‘ dBZ ) 


( 22 ) 


Substituting a dBZ of 25.5 into (22) results in (R ) =1.4 mm/hr. 

r min 

In summary, assuming radar parameters which are the same as the 
Seasat altimeter, a unity signal to noise ratio of -115 dBm results 
in a 25.5 beam averaged reflectivity factor which is equivalent to 
an estimated 1.4 mm/hr rain rate. 


6. DECONVOLUTION ON AN ARBITRARY REFLECTIVITY PROFILE 


In this Section we demonstrate the power of the deconvolution 
technique in recovering high resolution reflectivity profiles. We 
implement this technique employing deconvolution methods as de- 
scribed in Section 3 and radar data of a rain reflectivity profile 
acquired employing a ground based radar at Wallops Island, Virginia. 

In Figure la (solid curve) we show, as an example, a typical 
measured reflectivity profile over a distance interval of 90 km. 

This profile was obtained with a radar operating at S band (2.8 
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Distance traversed by sub-satellite point (km) 

Figure 1 (a) Measured reflectivity factor profile (solid curve) 

and beam averaged levels (dashed curve) . (b) Measured reflec- 

tivity factor profiles (solid line) and deconvolved profile 
(solid and dashed). (c) Gain squared beam pattern. 
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c- Hz), a Low elevation angle pointing (0.5° relative to the local 
horizon), and a 150 m range resolution (oulsewidth = 0.5 ps.'-c) 
(Goldhirsh, 1979). The profile corresponds to the rain day January 
2, 1979 ana was obtained at 9:34 a.m. local time. 


Shown in Figure l-(a) is a dashed curve which reoresents the beam 
averaged reflectivity factor with values samoled every 0.7 km. This 
is the reflectivity that would be measured assuming the footprint of 
the satellite borne radar (whose parameters are given in Section 4) 
sweeps by the reflectivity profile. This curve was calculated by 
injecting the true reflectivity profile into the integral (2), where 
a uniform reflectivity is assumed in the cross track direction. The 
profile was assumed periodic in the calculation of the beam averaged 
values in the vicinity of the 0 and 90 km distance intervals. We 
note that beam averaging considerably smooths the variations of the 
true profile. Ir fact, the smoothing appears so e-treme that it is 
difficult to visually correlate the beam averaged reflectivity with 
the actual profile. Also shown by the short dashed horizontal line, 
is the level of absolute reflectivity which would result in a unity 
signal to noise ratio. This corresponds to a beam averaged reflec- 
tivity factor level of approximately 25 dFZ. We rote the simulated 
measured beam average values exceed this lower threshold level . 


The right hand scale denotes the corresponding rain rate levels 
in mm/hr employing the empirical relationship given by , 


R = 4.0 X IQ’ 3 10 ( - 0S6 dBZ! 
r 


(23) 


This result was calculated from drop size distribution measurements 
mado at Wallops Island for the same rain event. We note that over 
the reflectivity interval, the rain rate levels are between .4 and 

7 mm/hr . 


In Figure 1- (c) (upper figure) is depicted the gain squared var- 
iation taken over the same abscissa scale as the actual reflectivity. 
It is interesting to note that considerable variation in the reflec- 
tivity profile exists over the 22.3 km footprint corresponding to 
the 6 dB down gain squared points. 

To implement the computer simulation, the Fourier and inverse 
transforms were replaced by their series approximations. Specifical- 
ly, the Fourier series of the gain-squared antenna pattern and the 
beam averaged reflectivity were computed; otherwise, the technique 
employed is the same as described in Section 4. 

In Figure 1- (b) the dashed line represents the high resolution 
reflectivity data that was recovered from the beam averaged reflec- 
tivity employing ‘•he deconvolution technique. Where no dashed line 
is observed, the recovered level coincides with the actual reflec- 
tivity values (solid line). We note that differences are evident 
mostly at the very low reflectivities, equivalent to negligible rain 
rates, and are due to computational imprecision. This computer sim- 
ulation confirms the thesis that the high resolution data can, in 
principle, be recovered from the beam averaged data. The above 
simulation was performed without injecting any contributions due to 
noise. The effects of noise on the resultant reflectivity profile are 
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examined in Section 7. 

It is interesting to note that the technique is able to re- 
cover reflectivity factors at levels below the threshold noise level 
for the beam averaged case (i.e., below 25 dBZ) . This may be attri- 
buted to the fact that the beam averaged measurement is a single mea- 
surement made by the receiver at a fixed point in space. On the 
other hand, the deconvolved measurement is one which employs a series 
of beam averaged measurements over ranges (in the example) , exceed- 
ing 100 km. The effective bandwidth of the system is considerably 
reduced since the effective integration time is substantially in- 
creased. 

7. EFFECT OF NOISE ON THE DECONVOLUTED SIGNAL 

In Figure 2 we demonstrate the effect of intrinsic noise on the 
deconvolved signal. In Figure 2- (a) is depicted the true reflec- 
tivity profile given by the solid line and the beam averaged echo 
(dashed line) . Superimposed on the beam averaged echo is the Pay- 
leigh fading noise introduced by the rain itself. As an example, we 
assume the capability of measuring 400 independent samples. The es- 
timate of the average power for such a case follows a Gaussian prob- 
ability density distribution (Goldhirsh, 19 n 'j) . A sample size of 
400 hence produces a standard error of .05* (reciprocal of the square 
root of the sample size) or approximately 0.2 dB rms statistical un- 
certainty. 

In Figure 2- (b) are shown deconvolved signals for two cases. 

The solid curve represents the deconvolved result derived from a 
beam averaged estimate corresponding to 400 independent samples as 
described in the above paragraph. As a basis of comparison, v/e 
show also the deconvolved result derived from the beam averaged es- 
timate employing an extreme sample size of 10^ independent samples. 

8. CONCLUSIONS 

We note that introducing noise to the technique seriously de- 
grades the deconvolved result. However, in spite of this degrada- 
tion, considerably more information is obtained regarding the magni- 
tude and scale dimensions of the rain employing the convolution 
method than is obtained by the beam averaged case. 

Although this method has been tailored for the measurement of 
rain employing a radar, its application to other types of satellite 
remote sensors should be explored. 
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